No spin-glass transition in the "mobile-bond" model 
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The recently introduced "mobile-bond" model for two-dimensional spin glasses is studied. The 
model is characterized by an annealing temperature Tq. On the basis of Monte Carlo simulations 
of small systems it has been claimed that this model exhibits a non-trivial spin-glass transition at 
finite temperature for small values of Tq. 

Here the model is studied by means of exact ground-state calculations of large systems up to 
A'' = 256^. The scaling of domain- wall energies is investigated as a function of the system size. 
For small values Tq < 0.95 the system behaves like a (gauge-transformed) ferromagnet having a 
small fraction of frustrated plaquettes. For Tq > 0.95 the system behaves like the standard two- 
dimensional ± J spin-glass, i.e. it does not exhibit a phase transition at T > 0. 

PACS numbers: 75.50.Lk, 05.70.Jk, 75.40.Mg, 77.80.Bh 



Spin glasses^ are the prototype model for disordered 
systems investigated extensively during the last three 
decades in statistical physics. These systems exhibit 
complex energy landscapes resulting in many interest- 
ing phenomena like glassy behavior and aging. Despite 
much effort, still many open questions exists. The ques- 
tion about the lower critical dimension of Ising spin 
glasses had been discussed for quite a whil o^i^i'^i^ . Now 
it is clear that in two dimensions non stable spin-glass 
phase at finite temperature exists^-^AMSi. This has mo- 
tivated the search for other two-dimensional spin-glass- 
like systems exhibiting a Tc > 0^- Recently D.K. Sunko 
has proposediSi a "mobile-bond" model where qucnched- 
disorder realizations of ± J spin glasses are created by an 
annealed simulation, allowing the bonds to move. The 
system is equilibrated at high temperature, followed by 
a quench to a temperature Tq. Sunko has performed 
Monte-Carlo simulations of systems up to size L = 16 
and claimed that for low quenching temperatures Tq the 
model exhibits a spin-glass transition at finite tempera- 
ture Tc > 0. 

In this rapid communication, the model is studied by 
means of exact ground-state calculations of large systems 
up to L = 256. The scaling of domain- wall energies^iii^ 
is studied as a function of the system size. It is shown 
here that the model exhibits no spin-glass transition at 
finite temperature. For small values of Tq < 0.95 the 
system exhibits ferromagnetic order, while at Tq a tran- 
sition to the normal two-dimensional spin-glass behavior 
is found, i.e. Tc = 0. 

The model consists oi N = Ising spins Si — ±1 on 
a square lattice with the Hamiltonian 

H = —''^ JijSiSj , (1) 

where the sum runs over all pairs of nearest neighbors 
and the ± J are quenched random variables. 

The realizations are prepared exactly in the same way 
as in Ref. 12. For each realization, first N bonds with 
strength +J and TV bonds with strength —J are dis- 
tributed randomly among all 2N bonds. Then the val- 



ues of all spins are set randomly to orientations Si = ±1. 
Next, an annealed Monte-Carlo (MC) simulation!^ is per- 
formed. This means, at each step either a spin is allowed 
to flip or two bonds incident to the same site are al- 
lowed to exchange their positions. Each choice occurs 
with probability 0.5. Each step is accepted with the usual 
Metropolis probability depending on the energy change 
according to the Hamiltonian JQ. First, the system is 
equilibrated at high temperature T = 5 for 1000 MC 
steps per spin (MCS)i^. Finally the system is quenched 
to T = Tq and simulated for further 1000 MCS. The re- 
sult is a realization of the disorder which can used for 
further treatment, here ground-state calculations are ap- 
plied. 

In greater than two dimensions, or in the presence of a 
magnetic field, the exact calculation of spin-glass ground 
states belongs to the class of NP-hard problemsiSiil. 
This means that only algorithms with exponentially in- 
creasing running time are known. However, for the spe- 
cial case of a planar system without magnetic field, e.g. 
a square lattice with periodic boundary conditions in at 
most one direction, there are efficient polynomial-time 
"matching" algorithms'^. The basic idea is to repre- 
sent each realization of the disorder by its frustrated 
plaquettesi^. Pairs of frustrated plaquettes are con- 
nected by paths in the lattice and the weight of a path 
is defined by the sum of the absolute values of the cou- 
pling constants which are crossed by the path. A ground 
state corresponds the set of paths with minimum total 
weight, such that each frustrated plaquette is connected 
to exactly one other frustrated plaquette. This is called 
a minimum- weight perfect matching. The bonds which 
are crossed by paths connecting the frustrated plaquettes 
are unsatisfied in the ground state, and all other bonds 
are satisfied. 

For the calculation of the minimum-weight per- 
fect matching, efficient polynomial-time algorithms are 
available22i2i. Recently, an implementation has been 
presented^^, where ground-state energies of large systems 
of size N < 1800^ were calculated. Here, an algorithm 
from the LEDA library— has been applied, which allows 
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a quick implementation. It was not necessary to go be- 
yond N = 256^ (with is much larger than N — 16^ in the 
original worliiSi) to obtain reliable results. 

To study whether an ordered phase is stable at fi- 
nite temperatures, the following procedure is usually 
appliediSiMiLiiiSiiS£. First a ground state of the system 
is calculated. Then the system is perturbed to intro- 
duce a domain wall and the new ground-state energy is 
evaluated. Typically, the system initially has periodic 
boundary conditions in both directions, and the pertur- 
bation involves replacing periodic by antiperiodic bound- 
ary conditions in one direction. The domain-wall energy 
AE is given by the difference of the two ground-state en- 
ergies. In case the model exhibits long-range ferromag- 
netic order at non-zero temperatures, the domain-wall 
energy, averaged over many independent samples, has to 
increase with system size. E.g. for a pure 2d ferromagnet, 
the domain wall consists of a straight line, resulting in 
AE ~ L. For a spin-glass, none of the ground states with 
periodic and antiperiodic boundary conditions has a pri- 
ori a lower energy. Hence one studies the absolute value 
of the domain-wall energy to detect whether the systems 
exhibits spin-glass ordering at finite temperatures. 

However we cannot apply the matching algorithm for 
boundary conditions which "wrap around" in both direc- 
tions. For this reason, here the periodic boundary condi- 
tions in the y-direction are broken for each realization (by 
setting the bonds connecting the first and the last row to 
zero). This has no influence on the fact of whether the 
systems orders or not because the change of the bound- 
ary conditions to create the domain walls occurs in the 
x-direction perpendicular to the open boundaries. 

Here system sizes L = 4, 6, 8, ... 192, 256 are consid- 
ered. For each size, 1000 independent realizations of 
the disorder were generated for quenching temperatures 
Tq = 0.1,0.9,0.95,1.0 and 1.5. Then ground states 
with periodic (P) and antiperiodic (AP) boundary con- 
ditions in x-direction were calculated using the exact 
matching algorithm, resulting in ground-state energies 
Ep resp. Ep^p. The change in the boundary conditions 
introduces a domain wall in each realization with energy 
AE = EIp ~ E^. 

First, we consider a very low quenching temperature 
Tq = 0.1, which was claimed in Ref. to exhibit a spin- 
glass transition at T^/J = 2.22(1) and no ferromagnetic 
order. In Fig.nthe distribition P{AE, L) over the disor- 
der of the domain- wall energies is displayed for different 
system sizes L. Clearly, the domain- wall energies grow 
strongly with system size, which is an indicator for fer- 
romagnetic order. Below (c.f. Fig. it is shown that 
indeed the disorder average (AE) grows linearly with L, 
as in the normal ferromagnet. But the model exhibits no 
global magnetic moment, as found already in Ref. 0. 
This is due to fact that 50% of all bonds are antiferro- 
magnetic. Nevertheless, the model behaves like a ferro- 
magnet. The reason is that the bonds are distributed in 
the system such that only few frustrated plaquettes are 
present. Hence, each realization can be mapped via a lo- 
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FIG. 1: Distribution P{AE,L) of domain-wall energies for 
Tq = 0.1 and system sizes L = 8, 32, 64, 128, 256. 



cal gauge transformation on a ferromagnet with a small 
number of antiferromagnetic bonds. This explains the 
fact that in Ref. |l^ the critical exponent of the correla- 
tion length found at Tc/J = 2.22 was indeed that of the 
pure ferromagnet. 




FIG. 2: Distribution of domain- wall energies for Tq 
system sizes L = 8, 256. 
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Next, a large quenching temperature Tq = 1.5 is con- 
sidered. In Fig. 13 again the distribution of domain-wall 
energies for different sizes are show. For large sizes, the 
distrubtions are centered around AE = indicating the 
absence of ferromagnetic order. Furthermore, the width 
of the distributions decreases slightly with increasing sys- 
tem size, which shows that spin-glass order is not stable 



3 



against thermal fluctuations. This is the usual situation 
found for the two-dimensional ±J spin glass^'® (having 
Tc/J^O). 
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FIG. 3: Fraction of frustrated plaquettes present after the 
quench to temperature Tq as a function of Tq. The line is a 
guide to the eyes only. 

To understand the behavior of the system better, 
next it is analysed as a function of Tq. The behav- 
ior is probably mainly determined by the fraction of 
frustrated plauqettes. In Fig. |31 the average fraction 
of frustrated plaquettes of the quenched realizations is 
shown as a function of Tq. This can be compared with 
the standard ±J random bond model (with parameter 
p (z [0, 1]), which has on average 2Np antiferromagnetic 
and 2N{1 — p) ferromagnetic bonds. This results in a 
average fraction 

X{ 4:x{l - x)[x^ + {1 - x)^] (2) 

of frustrated plaquettesS^. For the ±J model, a ferro- 
magnet spin-glass (Tc = 0) transition occur9£i2L2S422i2£ 
near p = 0.11 were a;/ (0.11) « 0.31. This corresponds to 
a quenching temperature Tq « 0.95, see Fig. |31 Hence, 
for a comparison, simulations near p — 0.11 for the 
±J model and near Tq 0.95 for the "mobile-bond" 
model have been performed. Furthermore both models 
were investigated for two other pairs of parameters ex- 
hibiting similar concentrations of frustrated plaquettes: 
Tq = 0.9; p = 0.1 and Tq = 1.0; p = 0.12. 

In Fig. 01 the mean value (AE) of the domain-wall 
energy is shown as a function of the system size for 
Tq = 0.1,0.9,0.95,1.0. For Tq = 0.1 a clear linear in- 
crease occurs, corresponding to a normal ferromagnet. 
For Tq — 0.9 the domain- wall energy still increases with 
system size. The resulting values are very similar to the 
domain- wall energies found at p = 0.1 for the ±J model. 



FIG. 4: Mean value of the domain-wall energy (A_E) as a 
function of system size for Tq — 0.1,0.9,0.95, 1.0 (displayed 
by lines). Also the domain- wall energy for the ±J random- 
bond model with concentrations p = 0.1, 0.11, 0.12 of the an- 
tiferromagnetic bonds is shown (symbols). 



For Tq — 0.95 and Tq = 1 the mean domain-wall ener- 
gies decrease as a function of the system size, hence no 
ferromagnetic order persists. In these cases the data is 
almost equal to the results for p = 0.11 resp. p = 0.12 
of the ± J model. Please note that it is not claimed here 
that e.g. Tq = 0.9 corresponds exactly to p = 0.1. But it 
seems certainly possible to chose p such that the results 
for both models agree exactly. 

In Fig.Othe corresponding results for the mean (| Ai?|) 
of the absolute value of the domain- wall energy is shown. 
For values Tq < 0.9 again an increase is observed, due to 
the increase of the mean (non absolute) (Ai?). If spin- 
glass ordering existed in a system, then (jAi?!) would 
increase with growing L, while (Ai?) has to decrease. For 
Tq = 0.95, 1.0 (|Ai?|) increases only for small system sizes 
(which may cause signs of a stable spin-glass phase when 
simulating only small systems), while it starts to decrease 
with L for larger values of L. Hence spin-glass order is 
destroyed for any finite temperature T > 0. Please note 
that again the results at Tq — 0.9,0.95,1.0 agree well 
with the results aX p = 0.10, 0.11, 0.12 for the ± J model. 

To conclude, in this work the recently proposed 
"mobile-bond" spin-glass model has been investigated. 
An exact ground-state matching algorithm has been ap- 
plied, allowing to study large system sizes like N = 256^. 
The model turns out to be mainly equivalent to the 
± J random-bond model, which has been studied exten- 
sively in the past. Hence, for low values of the annealing 
temperature Tq, the model (corresponding to small con- 
centrations p of the antiferromagnetic bonds in the ± J 
model) exhibits ferromagnetic order. The only difference 
is that the Sunko model exhibits no magnetic moment. 
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FIG. 5: Mean absolute value of the domain- wall energy 
(1AS|> as a function of system size for Tq = 0.1, 0.9, 0.95, 1.0 
(displayed by lines). Also the absolute value of the domain- 
wall energy for the ±J random-bond model with concentra- 
tions p — 0.1,0.11,0.12 of the antiferromagnetic bonds is 
shown (symbols). 



since by construction the number of ferromagnetic bonds 
equals the number of antiferromagnetic bonds. Both 
models can be mapped onto each other by local gauge- 
transformations, the characteristic parameter is the frac- 
tion of frustrated plaquettes. 



For larger values of Tq > 0.95 (corresponding to p > 
0.11) the model displays the standard behavior of a two- 
dimensional ±J spin-glass, hence no order for T > 
exists. To summarize, the "mobile-bond" model does 
not exhibit a finite-temperature spin-glass transition at 
any value of Tq, opposed to the claims made in Refii^i. 
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